A unifying model for several two-dimensional phase transitions 
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A relatively simple and physically transparent model based on quantum percolation and dephasing 
is employed to construct a global phase diagram which encodes and unifies the critical physics of 
the quantum Hall, "two-dimensional metal-insulator", classical percolation and, to some extent, 
superconductor-insulator transitions. Using real space renormalization group techniques, crossover 
functions between critical points are calculated. The critical behavior around each fixed point is 
analyzed and some experimentally relevant puzzles are addressed. 
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Two-dimensional phase transitions have been a focus of 
interest for many years, as they may be the paradigms of 
second order quantum phase transitions (QPTs). How- 
ever, in spite of the abundance of experimental and theo- 
retical information, there are still unresolved issues con- 
cerning their behavior at and near criticality, most simply 
exposed by the value of the critical exponent v (describ- 
ing the divergence of the correlation length at the transi- 
tion) . Below we discuss three examples which underscore 
these problems. First, in the integer quantum Hall (QH) 
effect, which is usually described within the single parti- 
cle framework, various numerical studies yielded a crit- 
ical exponent v ss 2.35 Q, in agreement with heuristic 
arguments [2] . Some experiments indeed reported values 
around 2.4 |3j, while other experiments reported expo- 
nents around 1.3 |4|, close to the classical percolation 
exponent v p = 4/3. Even more perplexing, some exper- 
iments claim that the width of the transition does not 
shrink to zero at zero temperature in contradiction 
with the concept of a QPT. Additionally, the observa- 
tion of a QH insulator |(| is inconsistent with the QPT 
scenario Consider secondly the superconductor (SC) 
- insulator transition (SIT), for which theoretical stud- 
ies suggest several scenarios. Similar to the QH situa- 
tion, some experiments yield a value v ~ 1.3 not 
far from the value v ~ 1, predicted by numerical sim- 
ulations within the random boson model, but closer to 
the classical percolation value. Other experiments, how- 
ever, yield v ~ 2.8 while some experiments claim an 
intermediate metallic phase |lfjj . As a third example, 
consider the recently claimed metal-insulator transition 
(MIT) |TH . The critical exponent is again close to 1.3 
[12t Il3| but the occurrence of such phase transition is in 
clear contrast with the scaling theory of localization. 

The basic question which naturally arises is then 
whether it is possible to unify these two dimensional 
phase transitions within a single theory and thereby re- 
solve some of the problems raised above. A hint toward 
an affirmative answer is gained by experimental indica- 
tions that percolation plays a key role in both the QH 



transition 0] and in the SIT 01- Moreover, its rele- 
vance to the MIT has been argued theoretically and ob- 
served experimentally [H El El The QH and the SIT 
have been treated within a percolation-like model in Ref. 

consisting of SC or QH droplets connected via quan- 
tum tunneling, similar in spirit to the model presented 
in [13. 

In the present work we address these issues using 
a physically transparent picture. A similar approach 
proved to be quite successful for developing a model en- 
coding the QH transition which can also describe some 
aspects of the SIT 18]. Introducing a new aspect, de- 
coherence, into the model, we are able to include the 
QH transition (or the SIT), classical percolation and the 
MIT all within the same phase diagram. Employing real- 
space renormalization group (RSRG) techniques, we first 
calculate the critical exponents. Most important, we em- 
ploy this phase diagram for understanding the possible 
critical features, and expose the physical conditions nec- 
essary to observe them. The crossover between different 
critical points, which may be explored experimentally is 
investigated and concrete experimental predictions are 
made. 

For the sake of self-consistence, the model used to de- 
scribe the QH effect is briefly explained here. In 
high magnetic fields, electrons follow equi-potential lines, 
which, at low Fermi energy, are trapped in potential val- 
leys. At zero temperature, transport through the system 
is due to quantum tunneling between such trapped tra- 
jectories. Within the model, this picture is mapped on a 
lattice such that each closed trajectory in a valley corre- 
sponds to a site, and quantum tunneling between (near- 
est neighbor) sites corresponds to a (quantum) link. As 
the Fermi energy rises and crosses the saddle point sepa- 
rating two isolated contours, these two trajectories merge 
into one. The transmission between these two sites be- 
comes perfect, and the link connecting them is consid- 
ered SC. The QH transition point occurs when such a 
trajectory spans the entire system, which, in the model 
corresponds to percolation of the SC links. Consequently, 



the QH problem maps onto a lattice, with random tun- 
neling amplitudes between its sites, and a finite fraction 
of perfect bonds, namely a mixture of quantum and SC 
links. The critical behavior of the QH is encoded by the 
scaling behavior near the percolation point of this lattice 
model. Calculation of the transmission through such a 
system, using two different approaches, yielded a diverg- 
ing localization length, with an exponent v « 2.4 ± 0.1, 
close to the results of other numerical estimates 0. 

The new element introduced into the model is dephas- 
ing. A fraction q of the links is attached to current con- 
serving phase disrupting reservoirs 0]. Physically, it is 
a consequence of finite temperature but can result from 
other mechanisms. An electron propagating along such a 
link enters the reservoir and loses its phase, before it, or 
another electron, emerges from the reservoir (there is no 
net current into the reservoir). Transport along this link 
is therefore incoherent. Transmission in the presence of 
such phase breaking resistors can be treated within the 
Landauer-Buttiker formalism. The model parameters are 
then p, the fraction of SC links (i.e quantum links whose 
transmission is unity), and q, the fraction of classical, 
or incoherent links. Evidently, the fraction of quantum 
tunneling links is 1 — p — q. The principal objective is 
then to construct a phase diagram in the {p, q) plane of 
parameters and to identify the various critical transitions 
in this plane. 

Our main tool will be the RSRG scheme on a square 
lattice. In order to gain experience and confidence we 
first perform calculations on the line q — which reduces 
to our previous, fully coherent model as a special case. 
Rescaling the lattice by a factor of 2, each 2x2 square 
is mapped onto a 1 x 1 square (Fig la). Given the trans- 
mission amplitude of each link in the 2x2 square, we 
evaluate the transmission through the whole unit, by at- 
taching leads to its left and right sides (see below). This 
will allow us to follow the distribution of the transmission 
amplitudes from one iteration to another, and eventu- 
ally evaluate the distribution for the full lattice. For the 
classical problem this calculation reduces to the calcula- 
tion of transport through a Wheatstone bridge (Fig. lb), 
where the symmetry of the square lattice has been used. 
For simplicity we use the Wheatstone bridge geometry for 
the quantum problem, and later verify that the addition 
of dangling bonds does not modify the critical behavior. 

The actual calculation is carried out using the scatter- 
ing matrix approach. Each link carries a left going and 
a right going channel. In accordance with the physics at 
strong magnetic fields there is no scattering in the junc- 
tions (valleys) and the edge state continues uninterrupted 
according to its chirality (see Fig. lc), while the scatter- 
ing occurs on the link (saddle point) itself. Each scatterer 
is characterized by its scattering matrix Si, namely a 
transmission probability T; and phases. For each realiza- 
tion of transmission probabilities on the larger lattice, we 




FIG. 1: The real-space renormalization group process on a 
square lattice, (a) The mapping of a 2 x 2 square onto a 1 x 1 
square, resulting in reducing the size of the lattice by a factor 
of 2. (b) The equivalent mapping for which the calculation 
has been done. For the classical problem this equivalence is 
exact, while for the quantum problem it neglects the effects of 
dangling bonds. Later in the calculation dangling bonds have 
been included with no observable change in the critical behav- 
ior, (c) The scattering matrix approach (rotated 90 degrees 
compared to (b)). At each junction each edge state propa- 
gates freely according to its chirality. Scattering is along the 
links (marked as partially filled rectangles across the links). 
It is easy to see that if there is a percolating path of perfect 
scatterers (scatterers with unit transmission), then there is a 
perfect transmission along the whole structure. 



can evaluate exactly the transmission through the whole 
lattice T({Ti}, (f>) , and thus follow its distribution from 
one iteration to another, allowing us to determine the 
fixed points and the RSRG flow. If one starts with a 
distribution of scattering matrices whose phases are uni- 
formly distributed, they remain so, and thus we will be 
only interested in the scaling of the distribution of the 
transmission probabilities, 

G W (T) = J ...J 5 (t - f(Ti, . . . ,T 5 ,<fi)^ x 
x G^-^iT^-'-G^-^iTn) dT 1 ---dT 5 d<j>, (1) 

where <p stands for all the independent phases. 
The initial distribution can be written as 

G {0} (T) = pS(T-l) + (1 -p)G {0} (T), (2) 

where p is determined by the fraction of saddle-point 
energies that are below the Fermi energy, p(ep) — 
J_ F F(e)de, with F(e) the distribution of the saddle 

point energies. &°^(T) is determined by the relation be- 
tween the transmission and the saddle-point energy and 
by F(e). The initial distribution G {0} (T) does not affect 
the critical behavior, as the distribution always flows to- 
ward one of three possible fixed-point distributions. 

Since the probability that the transmission is unity is 
determined by the classical percolation probability, one 
finds 

G^(T) = P n (p)S(T-l) + (l-P n (p))G™(T,p), (3) 
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where P n (p) is the classical percolation probability after 
n iterations, and G^ can be exactly expressed in terms 
of & n ~ 1 } . Clearly, if p < p c , P{p) flows toward zero, 
while for p > p c it flows to unity. At p — p c the distribu- 
tion G flows toward a two-peak fixed distribution, similar 
to previously obtained results for the QH transition [2(| • 
The critical behavior can be deduced by investigating the 
length dependence of the averaged transmission near the 
critical point and its collapse using T(p, L) = T[L/^(p)], 
where £(p) is an energy (or concentration) dependent lo- 
calization length. We find that diverges at p c = 1/2 
with an exponent v = 2.4±0.f, consistent with previous 
numerical calculations for the QH transition 0. 

Having demonstrated the power of the RSRG tech- 
nique for q = (no dephasing) we now treat the model 
at an arbitrary point (p, q) . Qualitatively, several phase 
transitions can be identified. For q = the QH transition 
is recovered as demonstrated above. For p + q = 1 (no 
quantum links) the point q—p c ~\ describes the classi- 
cal conductor-superconductor percolation transition. For 
the case of p = 0, all the links are either quantum (with 
T < 1) or classical. This is the model suggested in Refs. 
[l2l Hq | to describe the " apparent" metal-insulator tran- 
sition in two dimensions. The flow lines can be deter- 
mined without the full quantum calculations by noticing 
that the rescaled cell will be SC if a cluster of SC links 
percolate, while it will be an incoherent metal if there is 
percolation of classical links, and no percolation of SC 
ones. These conditions define the RSRG equations for 
the quantities p and q 1 
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2(f-p)V + 8(1- pfp 3 + 5 {l-p)p i +p 5 (4) 
q {I0p 4 - 20p 3 (l-q) + 2p (1 - q) [2 + 5 (1 - q) q] 
q(2 + 2q- bq 2 + 2q 3 ) + 2p (2 + 3q - 10q 2 + 5q 3 ) } 



superconductor 



QHT- quantum Hall transition 
PT- percolation transition 
MIT- metal-insulator transition 




MIT metal 
(q=0.5) 



FIG. 2: The phase diagram of the model in term of the pa- 
rameters q, the concentration of a perfect-transmission links 
and p, the probability of incoherent links. The phase dia- 
gram displays three different phase transitions between the 
three different phases. The broken line represent a possible 
experimental trajectory leading to a finite metallic phase. 

The full phase diagram thus contains three limiting 



phases, an Anderson insulator, a superconductor (or a 
QH phase) and an incoherent metal, and is depicted in 
Fig. 2. The fact that for p > i one gets a QH phase is 
straightforward, as once there is percolation of SC links 
there is an edge state propagating without backscattering 
through the system. The transition between an Ander- 
son insulator and an incoherent metal is more intriguing. 
Once p + q > \ an electron can propagate freely from 
one side of the system to the other, but necessarily goes 
through incoherent scattering, and thus this phase is an 
incoherent metal. For p + q < ^ however, there is no 
trajectory that allows the electron to traverse the sys- 
tem without quantum tunneling. These tunneling events, 
however, are intertwined with incoherent scattering, and 
thus the conductance of the system will be of the order of 
cxp(— L^/c;), where ~ q^ 1 / 2 is the average distance 
between incoherent scattering events. As q decreases, 
the conductance of the system goes to zero, ending up 
as an Anderson insulator. For finite p, however, the ex- 
istence of incoherent scattering is not expected to affect 
the quantization of the Hall conductance [2lJ ■ Thus one 
may characterize the phase with finite p and q, but with 
p + q < \ as a "Hall insulator". In agreement with pre- 
vious treatments of this phase, we also conclude that it 
only exists for finite dephasing (q > 0), as for q = the 
system flows to an Anderson insulator. Further investi- 
gation of this phase will be reported elsewhere. 

Beyond the characterization of the different phases, the 
model allows a direct calculation of the crossover between 
them. The line at p = \ allows, as a function of q, to 
study the crossover between the quantum SIT (or the 
QH transition) to the classical superconductor-conductor 
transition. Following the RSRG flow along this line, we 
find that the average transmission obeys 



T(q,p=±)=T o (l+0 



L 



(5) 



With L<p(q) defined above and (3 some nonuniversal con- 
stant. This function correctly reproduces the two fixed 
point behaviors, T = T at the QH transition and 
T ~ at the classical percolation critical point. The 

values of t, the resistance critical exponent, and v p ob- 
tained separately by RSRG are 1.33 and 1.42, respec- 
tively, giving t/v p ~ 0.93, in good agreement with the 
limiting behavior of the crossover function J^l. The ef- 
fect of incoherent scattering on the critical behavior may 
explain the different critical exponents observed in QH 
and in SITs, which mostly agree with one of the critical 
exponents associated with these two critical points. 

The QH to the MIT crossover has been demonstrated 
experimentally [2^ . The present formalism allows us to 
study this crossover theoretically, by following the con- 
ductance along the line p + q = \. The dependence of 
the transmission on length, for different values of q along 
this line is depicted in Fig. 3. There is a clear crossover 
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from a constant (QH behavior) to T ~ L t / Up (classical 
behavior), as before. The length scale determining this 
crossover (see inset of Fig. 3) is found to be the percola- 
tion correlation length, namely for L < £ p the system be- 
haves quantum mechanically, and classically for L > £ p . 
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FIG. 3: Crossover from the QH critical point to the MIT: 
transmission as a function of length (on a log-log scale) for 
different values of (p c — p) , demonstrating the crossover from 
the QH behavior (length-independent transmission) to a clas- 
sical power-law dependence. Inset: The crossover length £ co 
as a function of the (p c — p) on a log-log scale, with a slope 
v — 1.23, close to the percolation correlation exponent. 

All these observations have clear experimental rele- 
vance. As mentioned above, the QH to the MIT crossover 
has been investigated experimentally, but with no de- 
tailed investigation of the critical behavior. We thus pre- 
dict that the critical exponent along this line will be the 
classical percolation exponent, v p , and not the QH expo- 
nent, with two relevant length scales for any finite sys- 
tem. Similarly, one can experimentally investigate the 
quantum SIT to the classical superconductor - conduc- 
tor transition, for different temperatures. The theory 
predicts that once the dephasing length becomes smaller 
than the system size (i.e. q becomes nonzero), the critical 
behavior will correspond to classical percolation. More- 
over, we predict that the conductance at the critical point 
will vary from its universal value at the quantum transi- 
tion to a length dependence value according to Eq. (01 . 

Another result of the proposed phase diagram is the 
existence of an intermediate metallic regime in the QH 
and in the SITs. Imagine changing a physical parameter 
(e.g. density or magnetic field) along the broken line in 
Fig. 2. Then there is one transition from an insulator 
to a metal, and then from a metal to a QH liquid (or 
a superconductor), consistent with experimental obser- 
vations for the QH transition^] and the SIT [Io|. The 
critical behavior of both these transitions is determined 
by classical percolation. Thus the theory predicts that 
the QH critical behavior can only be observed in system 
where these transitions coalesce into a single transition. 
This can be tuned, for example, by lowering the tem- 
perature and changing the dephasing length. Thus the 



theory can explain the different exponents observed ex- 
perimentally and predicts a crossover between the differ- 
ent critical behaviors as a function of temperature. 

It should be noted that some aspects of this phase dia- 
gram are not dissimilar to the one presented in Ref . |23| . 
Here, however, we do not invoke any zero temperature 
dissipation. Rather, if, as expected, the dephasing length 
diverges at T — > (i.e. q — > 0) the theory predicts only 
two stable phases - the QH (or superconducting) phase 
and the Anderson insulator phase. 
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